CVM studies on the atomic ordering in complex perovskite alloys 
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Abstract 

The atomic ordering in complex perovskite alloys is investigated by the cluster 
variation method (CVM). For the l/3{lll}-type ordered structure, the order- 
disorder phase transition is the first order, and the order parameter of the 1:2 
complex perovskite reaches its maximum near x=0.25. For the l/2{lll}-type 
ordered structure, the ordering transition is the second order. Phase diagrams 
for both ordered structures are obtained. The order-disorder line obeys the 
linear law. 
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I. INTRODUCTION 



Complex perovskite alloys draw immense attention because of their exceptional dielectric 
and piezoelectric properties, which have great actual or potential uses.EH Examples are the 
so-called "super Q" mixed-metal perovskites with low lossesiS, such as Ba(Zni/ 3 Ta 2 / 3 )03 
(BZT) and Ba(Mg!/ 3 Nb2/3)03-BaZr03 (BMN-BZ), and the relaxor ferroelectrics (relaxors) 
with extraordinarily high values of the relative permittivity and the piezoelectric constantsii, 
such as Pb(Sc 1/2 Ta 1/2 )03-PbTi0 3 (PST-PT) and Pb(Mg 1/3 Nb 2/ 3)0 3 -PbTi0 3 (PMN-PT). 
An important aspect in complex perovskite alloys is the compositional atomic ordering, 
which is closely related to the desired properties. For example, it was indicated that the 
microwave- loss properties of 1:2 perovskite ceramics are very sensitive to the B-site cation 
ordering]! In relaxors, it is widely accepted that the presence of the nanoscale ordered mi- 
croregions is responsible for the relaxor behaviors. H Some investigations have been done 
on the atomic ordering of complex perovskites. Considering the electrostatic interactions 
between the nearest B-site ions, Zhang et al. presented an available quantitative criterion 
of the order-disorder transition and successfully explained many experimental results in 
A(B' 1 y 2 B' 1 ^ 2 )03.i Recently, Bellaiche and Vanderbilt used the electrostatic model to satis- 
factorily reproduce the various types of ordered structures in BMN-BZ and other complex 
perovskites.0 In addition, the first-principle calculations were conducted to describe the 
ordering in some cases.EUli 

It has been revealed00 that the Fourier transforms of the atomic interactions play 
important roles in the atomic ordering. In this paper, we adopt the parameters in the 
reciprocal space of atomic interactions and conduct a cluster variation method (CVM)@ 
calculation to investigate the order-disorder transitions in complex perovskites. 
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II. METHOD 



When considering the atomic ordering in complex perovskites, we can model (1 — 
£)A(B^B / 1 / _ c )C , 3-£AB" / 03 as a three-component system (1 — x)(B' c B'l_ c )-xB'" on a simple 
cubic Bravais lattice. c=l/3 for 1:2 perovskites such as BMN-BZ, and c=l/2 for 1:1 per- 
ovskites such as PST-PT. 

The driving mechanism responsible for the ordering is mainly the electrostatic 
inter act ion. i'0 Introduce c(R) = -1, 0, or +1 in the case of c = 1/2 if R is occupied 
by B', B'" or B". In the case of c = 1/3, introduce cx(R) = -2, 0, or +1 for B', B'" or B". The 
charge fluctuation of ions at R, Ag(R), is proportional to cr(R).0 The electrostatic energy 
between a pair of cations at R and R' is proportional to Ag(R)Ag(R') or cr(R)cr'(R / ), i.e., 

E elec (v, cx\ R, R') = W(R - R')a(R)a'(R / ), (1) 

where VF(R — R') is the effective-interaction energy. In a purely Coulomb pictured, 
W(R) oc 1/R. In the nearest-neighbor Coulomb picture!, VF(R) is equal to zero unless 
R is the nearest neighbor. Different from the two pictures, we do not make any assumption 
on the form of W(R) here. The energy of the system with a certain configuration {cx(R)} 
can be expressed as 

E (MR)}) = E + ~ £ W(R - R>(R)<r'(R'), (2) 

Z R,R' 

where E is a constant independent on {<j(R)}.The average value of the energy is (Eq is 
neglected) 

E= %(R)}»(R)}) 

{a(R)} 

= \ E ^(R-R / )a(R)a'(R')X^(R,R'), (3) 

where P({cx(R)}) is the statistic probability of the configuration {cr(R)}, and X aa /(R, R') 
is the pair occupation probability. 
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When the atomic interaction is long range, Eq. (3) is difficult to handle in CVM be- 
cause a large cluster should be used to contain the long-range interaction. Therefore, some 
approximation should be made. 

Expand X aa i (R, R') in a Fourier transform: 

X aa , (R, R) = Q*«' ( k > k ') ex P[*( k R + k ' R/ )] , (4) 

k,k' 

and substituting it into Eq. (3) yields 

N 

E = - £ V(k)aa'Q*AK-k), (5) 

(T,cr',k 

where V(k) is the Fourier transform of VF(R): 

V(k) =J2W(R)exp(ik-R). (6) 

R 

For a certain ordered structure, only a few Qcr CT '(k, — k) in Eq. (5) are nonzero under 
the single-particle approximation. 00 We assume that this property is still valid in CVM. 
The assumption is not absurd since in experiments there are only a few superlattices for a 
certain ordered structure. Thus Eq. (5) can be simplified. For example, for the 1/2(111}- 
type-ordered structure [Fig. 1(a)], all lattice sites are divided into two nonequivalent groups. 
Only the superlattice Q a(J i(k = 1/2{111}) are nonzero. Then the energy is reduced to 

iV 

E = -V(k = 1/2{111}) £ aa'Q^(k = 1/2{111}). (7) 

cr,a' 

By using the symmetry of the ordered structure, we can evaluate Q aa /(k = 1/2{111}) from 
a small cluster and deal with Eq. (7) in CVM. It is noted that V(k = 1/2(111}) involves 
any long distance interactions. For the l/3(lll}-type ordered structure, lattice sites are 
divided into three groups, 

E = f^i/aE™' [<3-'( k = 1/3{1H}) + Qaa' (k = 2/3(111})] , (8) 

— era' 

where V 1/3 = V(k = 1/3(111}) = V(k = 2/3(111}). 

Four-point basic clusters (Fig. 1) are chosen in our CVM calculation. 



III. RESULTS 



The l/2{lll}-type ordered structure occurs in both 1:2 and 1:1 complex perovskites. 
V1/2 = V(k = 1/2{111}) < is required since it is a necessary condition for the 1/2(111}- 
type ordered structure to appear. 

First, the ordering transition type is specified. The long-range order parameters (LRO) 
are defined as 

77=1^-^-1, (9) 

which represent the occupation-probability difference between nonequivalent sites, rj for B', 
B", and B w as functions of the reduced temperature T* are plotted in Fig. 2 when c = 1/3 
and x = 0.25. (The reduced temperature is defined as T* = k B T/\Vi/2\ in this figure and the 
following Fig. 3.) It shows that r] decreases gradually to zero when increasing temperature. 
So the ordering transition is the second order. For 1:1 complex perovskites (c=l/2), the 
transition is also the second order. 

Fig. 3 depicts the curves of r\ as functions of the composition x when c=l/3 and T* = 0.5. 
It shows that t/b' increases when decreasing x, while rfe" and T]b"' reach their maxima near 
x = 0.25. The appearance of the maximum in LRO is consistent with the experimental 
findings in BZT-BZ0 and previous theoretical workH3. For 1:1 complex perovskites (c = 
1/2), rjB' and 7/b» are equal to each other and increase with decreasing x, but T]b'" keep equal 
to zero in all range of x. 

The critical temperatures of the phase transition, T c , can be determined at different x 
values. Approximately, T c decrease linearly with increasing x. 

The l/3(lll}-type ordered structure is observed only in 1:2 complex perovskites. The 
calculated curves of the long-range order parameters rj as functions of the reduced temper- 
ature T* when £=0.25 are shown in Fig. 4. (In this figure and thereafter, T* is defined in 
T* = ksT /\Vi/s\ where V1/3 < 0.) It can be seen that LRO of the l/3(lll}-type ordered 
structure suddenly drop to zero at the transition temperature, which indicates that the 
transition is the first order. 

5 



The 1/2{111}- and l/3{lll}-type ordered structures are separately discussed above. In 
actual cases, which ordered phase to appear is determined by the free energy, i.e., at a certain 
composition, the ordered phase with the minimal free energy will become the equilibrium 
phase and can be observed in experiments. 

In the investigations above, the model is parameter-free in the sense that V1/2 and V1/3 
define the temperature scale. When we consider the 1/2{111}- and l/3{lll}-type ordered 
phases at the same time, a ratio V1/2/V1/3 should be introduced. For the purely Coulom- 
bic interactions^, \V\/2\ ~ IV1/3I; for the nearest neighbor interaction!, [ V1/2 j = 2| V1/3 1 . 
Adopting I Vx/2 1 = 1.4|Vi/3|, we plot in Fig. 5 the free energy curves of the 1/2{111}- and 
1/3(111 }-type ordered structures for c = 1/3. It shows that the 1/3(111 }-type ordered phase 
appears at small x values, and the 1/2(1 ll}-type phase favors larger x values, which agrees 
with results in BMN-BZ0. The boundary of x values for the 1/2(111}- and l/3(lll}-type 
ordered phases is determined by the crossing of the curves. The combined phase diagram of 
the system is shown in Fig. 6. The phase region of the l/3(lll}-type ordered phase near 
x = is located inside the region of the l/2(lll}-type ordered phase. 

When I V1/3 1 decreases, the phase region of the 1/3(1 ll}-type ordered structure gets 
smaller and smaller, and vanishes at last. This may be the cases in PMNiffl and PMTH 
in which only the l/2(lll}-type ordered phase is observed. Of course, the source of the 
small IV1/3I cannot be solved in the theoretical framework of this paper. The effective energy 
parameters are affected by the A-site and O-site ions. Perhaps the different effects of the 
Pb and Bail can provide some clues on this problem. 

For 1:1 complex perovskites (c = 1/2), the chemical formula is favorable to form the 
l/2(lll}-type ordered structure, so the situation is different from that of c = 1/3. The 
free energy curves are depicted in the inserted graphics of Fig. 9 when c = 1/2. The curve 
of the l/3(lll}-type ordered structure wholly lies above that of the 1/2(1 11 }-type phase. 
It suggests that the l/3(lll}-type ordered phase is impossible (or difficult) to form in 1:1 
complex perovskites such as Pb(Sci/ 2 Tai/2)03i and Pb(Mgi/ 2 Wi/2)03il. 
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IV. CONCLUSIONS 



The 1/2{111}- and l/3{lll}-type atomic ordering in complex perovskites have been 
studied by the cluster variation method. For the l/2{lll}-type ordered phase, the tran- 
sition is the second order, and the order parameter reaches its maximum near x=0.25 for 
1:2 complex perovskites. For the l/3{lll}-type ordered phase, the transition is the first 
order. The order-disorder transition line is found to obey the linear law. The phase diagram 
containing both ordered structures is obtained. 
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FIGURES 

FIG. 1. Lattice site occupations and basic clusters used in CVM for: (a)the l/2{lll}-type 
ordered structure; (b)the l/3{lll}-type ordered structure. 

FIG. 2. LRO of the l/2{lll}-type ordered structure as functions of reduced temperature 
T* = k B T/\V 1/2 \ at x = 0.25 and c = 1/3. 

FIG. 3. LRO of the l/2{lll}-type ordered structure as functions of x at T* = 0.5 and c = 1/3. 

FIG. 4. LRO of the l/3{lll}-type ordered structure as functions of reduced temperature at 
x = 0.25 and c = 1/3. The reduced temperature is defined as T* = k B T/\Vi/ 3 \. 

FIG. 5. Reduced free energies of the l/2{lll}-type (solid line) and l/3{lll}-type (dashed line) 
ordered structures as functions of x when c = 1/3. V 1 / 2 = —1-4, Vi/ 3 = —1.0, and k B T = 0.6|Vi/ 3 |. 
Inserted graphics is the case of c = 1/2. 

FIG. 6. Combined phase diagram when c = 1/3. a, (3, and 7 designate disordered, 
l/2{lll}-type and l/3{lll}-type ordered phases, respectively. V 1 / 2 = —1-4, Vi/ 3 = —1.0. The 
reduced temperature is defined as T* = k B T /\Vi/^\. 
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